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ABSTRACT 

Explaining variability observed in the accretion flows of black hole X-ray binary systems remains 
challenging, especially concerning timescales less than, or comparable to, the viscous timescale but 
much larger than the inner orbital period despite decades of research identifying numerous relevant 
physical mechanisms. We take a simplified but broad approach to study several mechanisms likely rel- 
evant to patterns of variability observed in the persistently high-soft Roche-lobe overflow system LMC 
X-3. Based on simple estimates and upper bounds, we find that physics beyond varying disk/corona 
bifurcation at the disk edge, Compton-heated winds, modulation of total supply rate via irradiation 
of the companion, and the likely extent of the partial hydrogen ionization instability is needed to 
explain the degree, and especially the pattern, of variability in LMC X-3 largely due to viscous damp- 
ening. We then show how evaporation-condensation may resolve or compound the problem given the 
uncertainties associated with this complex mechanism and our current implementation. We briefly 
mention our plans to resolve the question, refine and extend our model, and alternatives we have not 
yet explored. 



1. INTRODUCTION 

The X-ray spectrum of black hole X-ray binaries 
(BHXRBs) often shows variability in intensity and hard- 
ness on timescales of order the viscous timescale, but 
much larger than the innermost orbital period, including 
the well-known "q" -dia gram hysteresis patterns traced 



by transient BHXRBs (Fender et al. 1999[ Hgman fc 
Belloni| |2005| |Done et al.||2007l and seFEgjl]). The 



properties ot such variabihty may also evolve over mul- 
tiple viscous-timescale cycles. The high-soft (but sub- 
Eddington) and quiescent intensity-hardness l imits are 
understood as manifestations of the thin-disk dShakura 



&: Sunyaev 1973 1 and radiatively -inefficient adyection- 
dominated accretion flow (ADAF; Narayan fc Yi| [T995al 
accretion limits while states between are typically under- 
stood as some evol ving combination of disk and ADAF 



flows (Chakrabarti fc Titarchuk 1995 
Nandi et al 



_^ Esin et al.||1997 

'I'heoretical work has begun in this 



ii2012| ). 

regime, but still cannot fully explain observations, espe- 
cially regarding viscous timescales where detailed simu- 
lations require prohibitively many time steps and steady- 
state assumptions lose validity. This motivates us to de- 
velop theory in more detail starting with a system like 
LMC X-3, whose behavior is more constrained than that 
of the transients, but still exhibits substantial variability 
that is quantitatively challenging to explain. 

In the current paradigm, transient BHXRBs are those 
systems where the outer disk reaches temperatures low 
enough to trigger the partial hydrogen ionization insta- 
bility (PHII) in which accretion proceeds via cycles as 
viscosity alternately concentrates m ass i nto rings and d if- 
fuses it inward (see Cannizzo||1998 and LasotapOOl for 
a review). The part ot the cycle where viscosity concen- 
trates mass leads to the quiescent phase where any emis- 



sion is presumably powered by some fraction of the flow 
that escapes the viscosity trap. Eventually, density in- 
creases enough to raise disk temperature above the tran- 
sition point in viscosity, and the sudden jump in accretion 
rate leads to a rise in luminosity (often several orders of 
magnitude) while still in the hard state, followed by a 
softening of the spectrum at this same, peak luminos- 
ity (the vertical and horizontal shifts by the dashed line 
in fig [I]). Jet emission is seen to s hut off as the system 
enters the extreme high-soft state (Fender et al. [1999 ). 
Afterward, the system typically makes partial transitions 
in hardness (with small, erratic shifts in luminosity) on 
sub-viscous timescales and associated with intermittent 
jet emission on top of a secular decline in luminosity. 
Eventually the system transitions completely back to the 
hard state and then finishes fading back into quiescence. 

There are a handful of persistent BHXRBs and black- 
hole-candidate X-ray binaries that do not execute such 
extreme quiescence-flaring cycles, and consistent with 
the transient paradigm they appear to avoid the PHII 
through an appropriate combination of disk size, l umi 
nosity (for disk irradiation) , and mass supply rate ( Co 



riat et a l.||2012[). Despite the label, these systems may 
still show significant variability and a variety of behav- 
iors. Cygnus X-1 appears to slide between canonical 
high-soft and l ow-hard states wit h some scatter but no 
clear hysteresis ( [Smith et al.|2002[ hereafter SHS02). The 
black hole candidates GRS 1758-258 and IE 1740.7-2942 
exhibit transient-like hysteresis but seldom reach the soft 
or quiescent limits (SHS02). LMC X-3 is typically bright 
in soft X-rays, but shows some hysteresis that circulates 
in the opposite direction to transients (see fig[ll, and 
does occasionally transition completely t o the hard state 
( [Wilms et all[200T| |Smale fc Boyd|[2012l ) . 
Accounting tor irradiation ot the outer disk, LMC X-3 
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Fig. 1. — Shaded points show a typical spectral evolution cycle of 
LMC X-3 (first episode in fig[5l with darker points corresponding to 
earlier observations, showingthat it winds around in the opposite 
sense of typical transients' cycles shown schematically in the gray, 
dashed "q" or "turtle-head". 



accretes at rates high enoug h to usuaUy avoid th e ioniza- 

though 



tion instabihty completely (Coriat et al. 2012 



the outer disk likely becomes susceptible tor the deeper 
drops in disk luminosity, and almost certainly for rare, 
complete state transitions. Such high accretion rates are 
explained via Roche-lobe overflow (RLO); optical mea- 
surements of the companion's spectral ty pe that take disk 
irradiati on into account indica te B5IV ( |Soria et al.|2001 1 
or B5V ( Val-Baker et al.|2007 ) spectral type. Such a star 
will fill the Roche lobe tor a 1.7 day orbit around a black 
hole with mass equal to the recent 9.5 Mp, lower bound 



(|Val-Baker et al."2007) . Furthermore, |Soria et al | ( |2001 1 
have argued that feeding the observed \x.-ray luminosities 
via winds would lead to column densities far higher than 
measured. Although its distance precludes direct mea- 
surement or exclusion of radio jets, the jet quenching ob- 
served in transient BHXRBs as they approach the high- 
soft state also suggests that LM C X-3 does not usually 
possess strong jets (Fe nder et al.| 1998). Thus, if LMC X- 
3 shares similar variability mechanisms with transients, 
while being far less prone to the ionization instability and 
jet outflows, then LMC X-3 provides a more controlled 
setting to study such mechanisms. Below we list the ma- 
jor mechanisms we have considered so far in modeling 
LMC X-3, also summarized in fig[2l 

For RLO systems like LMC X-37nard X-rays from the 
inner accretion disk can lead to supply-rate modulation 
(SRM) from the companion by inflating the companion's 
atmosphere to increase the density of gas at the LI La- 
grange point, as well as the a,rea and pressure behind the 
nozz' 



grange point, as well as tne a,rea ana pressure Denmd tne 
nozzle ([Lubow fc Shu||1975[ Meyer & Meyer-Hofmeister 
1983| ). We elaborate in ^3.1[ 

Xs the stream of gas from the companion dissipates 
energy and spirals onto the circularization radius, it may 
encounter the edge of the viscously spreading disk and 



from this point the flow can undergo disk corona bifurca- 
tion (DCB) as some fraction can efficiently shock, cool, 
and join the disk while some fraction may stream past 
the thin disk edge and maintain its vir ial temperature 



Hessman 1999 



Armitage & Livio'1998). 



W arping of the outer disk , whether driven by ir radia- 



tion ( Pringle 1992; Ogilvie & Dubus 2001; Foulkcs et al. 



2010| ) or a fift force ( |Montgomery fc Martin„2010p can 
attect the accretion flow by changing the d ensit y profile 



that the disk presents to the RLO stream (^3.2|, and by 
varying how the companion is exposed to or shadowed 
from inner disk X-ra ys. For LMC X-3 spe cifically, |Ogilvie| 



fc Dubus|p00T| ) and |Foulkes et aL| ( |2010[ ) predict that the 
system is potentially unstable to irradiation-driven warp- 
ing. Because the dominant long-term effects of warping 
are through bifurcation, and because we will lump them 
together as a boundary condition in our simulations, we 
combine discussion of them into §3.2| 

X-rays from the inner disk can Compton heat gas in the 
disk atmosphere and any corona at large radii above the 
local virial temperature thus dri ving a Compton-heated 
wind (CHW), discussed in ^3.3 which is not only im- 



portant for removing gas, but for removing hot corona 
that might otherwise help evaporate the disk or condense 
further inward (item (EC) below). As noted, for large 
enough disks and insufficient irradiation, a finite strip in 
the outer disk becomes susceptible to the PHIL For now, 
we do not treat it in any detail, but discuss th e lik ely 
extent and manner of its effects in LMC X-3 in 33. 4| 

If the disk and corona are coupled thermally, then the 
disk and corona may exchan ge mass through e vapora - 



tion and condensation (EC), 
hereafter MP07) provide a t 
numerical treatment, and Liu et al 



Mayer fc Pringle] (( |2007[ ) 
lorou g h introdu ction and 



j Liu et al.[ hereafter 



LTMHM07) and |Meyer-Hofmcister et al.| ( |2009[ ) discuss 
more applications and provide the steady-state prescrip- 
tion for our modified method. Through EC, extant disks 
will tend to preserve the soft state down to lower lumi- 
nosities via Compton-cooling-driven condensation, pro- 
viding a natural explanation of why BHXRBs return to 
the hard state at lower lumin osity and thus show hystere- 
sis ( Meyer-Hofmeister et al.|i2009, focus on this asp ect). 
We discuss other interesting effects possible in p.5| 

We will restrict our focus to an alpha-viscosity pre- 
scription for the disk. For the present work describ- 
ing long-timescale variability over accretion rates typical 
to LMC X-3 where uncertainties regarding conventional 
mechanisms still loom large, we consider this perfectly 
adequate, but acknowledge the possibility of more in- 
trinsic variability mechanisms (§5]). 

In S|2] we review key features of LMC X-3's accretion 
behavior, summarize how we infer the innermost disk 
and corona/ADAF accretion rates from the X-ray data 
(additional details are provided in the Appendix), and 
critically exa mine the quafi tatively simple bifurcation- 
only model in 'Smi th et^ (?, hereafter SDS07). The 
latter motivates ^ in which we furnish additional detail 
on the mechanisms as listed above, including estimated 
constraints on their effects and thei r cu rrent level of im- 
plementation in our modeling. In j ]4.1| we argue that a 
model including mechanisms besides EC cannot explain 
the data, but does best when given an unreasonably small 
disk radius and very large variations in c orona-disk ratio 
at this boundary. We then show in ^4.2 how EC may ef- 



Fig. 2. — A cartoon of the primary mechanisms cons idered for driving long timcscale changes observed in i nner dislt and corona accretion 
rates. These are: supply-rate modulation (SRM, i ]3.1| from the companion, disk-corona bifu rcat ion (DCB, j ]3.2| l at the disk's edge , wh ere 
disk warping may affect SRM and DCB, the partial hyd rogen ionization instability (PHII, j ]3.4| , Compton-heated winds (CHW, ^3.3^ at 
large radii, and evaporation and condensation (EC, i]3.5|l exchanging mass between disk and corona. 



fectively recreate such seemingly ad-hoc conditions, but 
how it may also imply behavior inconsistent with obser- 
vations, including an extremely easily triggered "sympa- 
thetic" mode where the innermost disk and corona accre- 
tion rates rise and fall simultaneously. We briefly review 
the results and caveats of the current model, and state 
our current plans to resolve the question in ^|5] 

2. LMC X-3 AS PROTOTYPE 

Besides simplifying initial modeling as discussed above, 
LMC X-3 offers additional practical advantages. The 
Rossi X-ray Timing Explorer {RXTE) monitored LMC 
X-3 for over 16 years, and at least five of those in- 
clude observations each about a kilosecond long taken 
roughly twice a week, thus providing a long, uninter- 
rupted history of accretion with sufhcient resolution 
at the timescales we seek to study. Also, the X-ray 
blackbody component, when present, tracks the Stefan- 
Boltzmann law fairly well (see fig|3]) indicating that inner 
disk geometry (i.e. truncation, warping) changes fairly 
little, and that the corona optical depth, Tc, is small, sim- 
plifying estimates of the inner corona accretion rate, SJl^. 

Unless specified otherwise, we will use symbols SJt^ (9^c) 
as shorthand for disk (corona) accretion rates at the inner 
disk radius. Rid, and reserve the italic-face for general, 
local accretion rates = Mj^{R,t), = M^{R,t). 
Also, unless otherwise indicated, we will use the follow- 
ing system parameters: black hole mass, M^^=1QMq, 
companion mass M*=5M(^, orbital period Pgyg=1.705d, 



inclination and system distance, dj,^=48 kpc 



van 



der Khs et al. f985 Val-Baker et al. 20071, which also 



imply a circularization radius of -R^irc ~ 2-7 x lO-'^^cm. 

We first fit individual RXTE spectra with a disk black- 
body and a power law of fixed photon index, Fpii = 2.34 
(as in SDS07) with total ab sorption o f fixed column den- 



sity riH = 3.8 X fO-^^cm"^ (Page et al. 2003 ), using th e 



wabs*simpl*diskbb models in XSPEC ( jArnaud 1996). 
To systematically identify transitions to the low-hard 
state we looked for cases where the first fitting gave 
reduced > 1.1, and refit these with a wabs*(plaw) 
model where the power-law index is not frozen. Reas- 
suringly, spectra identified this way were fit better with 
fewer parameters, and are also typically preceded by ob- 
vious declines in the blackbody component (fig|4]). 

For low Tc one can describe the flows qualitatively by 
taking 9)l^(i) ~ T^^ and 9)tc(^) proportional to the ra- 
tio of power-law to blackbody count fluxes (as in SDS07, 
though there the disk central temperature was confused 



Md hlOVs^O.lMEdd] 




Fig. 3. — Fitted blackbody-component fluxes plotted against fit- 
ted temperature with a pure T'* curve in faint gray for comparison. 

with the effective temperature giving ^^i^) ^ ^bb^^)- 
We obtain absolute normalization for 971^ by fixing Rid 
and comparing observed and predicted fluxes in the high 
state where agreement should be best, while for 9Jt^, we 
obtain an estimate based on the simple Tc and a typical 
ADAF solution, and check this against a more detailed 
calculation. We relegate the details to the Appendix 
to focus on a general description of accretion behavior 
(flg§. 

From figlsl one can see that the SJl^ "turns on" in 
pulses (referring to the secular month-long features and 
not the jagged week-long sub-pulses) roughly a vis- 
cous timescale apart and slightly shorter in duration, 
and that these pulses tend to antic ipate drops in 93T^ 



This trend was already noted in ( jSmith et al. 20071 



based on inferred qualitative accretion rates, and led 
the authors to posit a "bifurcation-only" model where 
a fairly-constant total supply rate (M^) is split far from 
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Fig. 4. — Fitted power-law (top panel) and blackbody (bottom panel) components of LMC X-3's X-ray flux since 53436.1 MJD. Diamond 
points in the top panel mark observations categorized as the pure hard state by the criteria in |2] 
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Fig. 5. — Inferred accretion history since 53436.1 MJD, barring times LMC X-3 was observed in the low-hard state (gray vertical lines) 
where inferring accretion rates is more ambiguous. Note the vertical labels refer to inner accretion rates here. In the top panel, solid trace 
shows simple, direct estimate of 971^, as well as results (points) of a more detailed method and arrowheads indicate points where the detailed 
method required abnormally high OT^, (see jjA|. Overall, one can see trend for OT^, to pulse "on" quasi-periodically and anticipate episodic 
drops in 931^ 



the black hole between non-interacting quickly-draining- 
corona and slowly-draining-disk components. Our nor- 
malization estimates for 971^(0 suggest that for any given 
episode there is generally insufficient total mass in a 
pulse to explain the associated EDt^ drop. Even if 
our overall normalization is off, we still found that scaling 
^ffl^ to conserve mass for one episode does not work very 
well for other episodes. This mass-conservation problem 
motivated considering mechanisms that can adjust the 
total supply rate, remove mass, and/or exchange it be- 
tween disk and corona flows. 

The secular evolution of the episodes on super- viscous 
timescales also lends itself to interpretation as multi- 
ple mechanisms acting on similar timescales effectively 



generating "beat-frequency" behaviors. The simple al- 
ternative of some mechanism(s) acting on super-viscous 
timescales coupled to viscous-timescale variability mech- 
anisms lacks good candidates for the former. Nuclear 
evolution is too slow, we do not expect significant mag- 
netic cycles from a companion with a radiative outer en- 
velope (but keep the possibility in mind regarding other 
systems), and the inferred mass ratio in LMC X-3 is too 
high for slowly-gro wing tidal resonances to be significant 
( Frank et al. [2002 ) . Furthermore, based on the observed 
1, thi 



inclination, the warps would have to reach heights of 30° 
relative to the orbital plane, and survive the severe drops 
in to exhibit precession effects if irradiation-driven, 
which poses difficulties if LMC X-3 is o nly marginally un- 
stable to irradiation-driven warping as |Ogilvie fc Dubus"] 



5 



( |2001[ ) suggest. 

'i'he disk component in LMC X-3 tends to fall and 
recover more rapidly for larger drops than for shallow 
drops, a trend quanti fied in SDS07 and rec ently over an 
expanded data set in |Smale fc B oyd (2012). This aspect 
is qualitatively consistent with a bifurcation-only model 
given sufficient variation in the amplitude and duration 
of a drop at the outer edge - sensitivity to duration for 
a single input amplitude can be seen in figures 7&8 of 
?. However, using their analytical machinery, with and 
without crude representation of SRM and outfiow effects, 
we will later show that rough quantitative agreement 
with observations of LMC X-3 requires inputs that are 
extremely unlikely without additional physics ([ ]4.1[ ). An 
interesting exception to the usual of DJl^ pulses heralding 
steep Wl^ drops is the small drop in at llOOd into 
fig ji] not associated with any 971^ pulse above the typical 
noise level. 

Inferring accretion rates in the absence of the disk 
component introduces additional parameters and uncer- 
tainties, but we wish to make a few relevant observa- 
tions while we work on a more definitive analysis of the 
hard state. The disk component drops and recovers on 
timescales of days in transitions into and out of the hard 
state, and tends to return more quickly than it decays 
when LMC X-3 is at its "hardest" in our data, circa 
the 1500d mark in bottom of figj5| consistent with the 
notion of an extant inner disk preserving itself through 
condensation. Also, the power-law component tends to 
increase before failed and successful disk restarts, which 
may physically correspond to the inner edge of a trun- 
cated disk moving inward to provide more and hotter 
seed photons, and/or rapid condensation. 

3. VARIABILITY MECHANISMS CONSIDERED 

Though the basic physics of companion irradiation and 
streaming are simple, the dynamics are potentially com- 
plicated to initialize and implement in detail, especially if 
the outer disk warps. However, we can estimate bounds 
on both mechanisms individually, and because they sit 
at the edge of the accretion fiows, we can lump them 
into a manual boundary condition for now and still de- 
rive meaningful results. Compton-heated winds can be 
launched a bit further inward, but can be described fairly 
well by simple analytical functions of radius and X-ray 
luminosity assuming that the corona is easily replenished, 
and thus we can quickly obtain upper bounds on CHW 
losses. 

Evaporation-condensation can depend sensitively on 
disk and corona conditions at all radii making it the least 
amenable to simple estimates, and as noted earlier this 
same strong dependence on the system's state can natu- 
rally engender hysteresis. EC also allows the disk com- 
ponent to vary more substantially and more rapidly by 
evaporating disk material interior to the circularization 
radius, but this evaporated disk material can also con- 
dense further inward much faster than inner disk condi- 
tions change, potentially to the point that 971^ rises and 
falls simultaneously with 9jl^. This "sympathetic" accre- 
tion mode can be seen in the more detailed simulations 
of Mayer and Pringle (their fig . 8) and in many cases we 
simulated (e.g. figures 8|11|12 1, hut is effectively absent 
(or negligible) in our observations of LMC X-3, and thus 



primarily poses a challenge to our basic EC model. 

3.1. SRM Estimates and Remarks 
The total supply rate of mass through the LI nozzle 
Mg, will scale with the product of local gas density Pl^, 
speed at which gas streams through the nozzle (roughly 
the local sound speed c^), and area of the nozzle A„ 
where the latter has width and height roughly equal 
to the isotherm al scale height in th e local tidal field, 
Hl^ sa c^jnl,^ ( |Lubow fc Shul|1975| ). Under X-ray ir- 
radiation, each layer ot the atmosphere will tend to heat 
up until it emits the intrinsic stellar flux plus the incident 
X-ray fiux at that altitude. Due to the very steep tran- 
sition in density at the photosphere, we find most of the 
X-ray energy is deposited in a thin layer there, which 
we will take to be infinitesimally thin for now. Thus, 
the modulation with respect to a given reference state 
as a function of stellar temperature T^, effective incident 
X-ray luminosity L^^cS^ gravity-darkened stellar luminos- 
ity offj and distance bet ween LI and the photosphere 
dLi — dph is given by (e.g. Meyer & Meyer-Hofmeister 



19831 



where 



yrcf 



3/2 



exp 



dhi 



rrrcf 



T] 



xef 



rcf 



1/4 



T rcf 



(1) 



(2) 



Short of solving the structure of the stellar envelope in 
the Roche-lobe potential under time-varying irradiation, 
we can estimate the extent of SRM by computing the ra- 
tio of effective incident-to-intrinsic luminosity. One can 
make a simple estimate by computing the effective grav- 
ity at a point sitting about halfway between the nozzle 
and the pole of the companion giving L„ off/i* ~ 0.68, 
and also use the inclination of this point relative to the 
inner disk to get the fraction of X-rays emitted into this 
latitude, cosf3x = 0.28, yielding 



max(La,_eff) < 0.3Lgjj x cos /J^r x 



47ra2 



< IO^Lq X 0.28 X 0.02 w lOOLg 



(3) 



The companion's effective stellar luminosity falls within 
^500-IOOOLq based on the reported bolometric stellar 



luminosity 800-1600L(;) (Soria et al. 2001). More care 



fully integrating the incident-to-mtrmsic ratio over the 
irradiated face (again, with gravity darkening) agrees 
closely with this simple estimate as the projected area 
and fraction of disk flux fall concurrently with (and faster 
than) the effective gravity toward LI. 

For irradiation operating alone, choosing the maximum 
observed luminosity as the reference point in eqnjljwould 
permit drops to ~50% of the observed maximum and 
only if the X-ray source were turned off completely, but 
this estimate is still fairly sensitive to companion temper- 
ature. Harder and more isotropic X-ray flux from a hot 
corona may enhance modulation, but for LMC X-3 the 
maximum observed power-law flux is barely a flfth that 



6 



of the disk, roughly equal to the projection factor reduc- 
ing inner-disk flux onto the companion. However, even 
if much deeper drops are possible, and irradiation-driven 
warping or some other mechanism were included to pre- 
vent the system from settling into a permanent steady 
high-soft state, the fact that SRM affects the flow at the 
very boundary means that any ch ange s it introduces will 
suffer severe viscous dampening (j j4.f [ ). Altogether, this 
suggests that SRM is significant, but certainly cannot 
explain the steep 971^ declines by itself. 

Furthermore, we consider this simple model's predic- 
tions of the SRM magnitude an upper bound in light of 
as detailed two-d imensional hydrodynamic s imulations of 
the envelope by Viallet & Hameury (2007). They find 
that irradiation will still drive gas toward the nozzle, 
but the gas will also have ample time to cool down as 
it crosses the disk's shadow. They note that because 
they do not solve for perpendicular velocity it may ex- 
ceed their estimates near the nozzle, and we remark that 
warping of the outer disk might reveal more of the com- 
panion's equator and nozzle and negate the effects of 
cooling. For our disk/corona simulations, we ignored the 
delay between irradiation and changes in since we es- 
timated the sound-crossing time of the envelope near LI 
to be wIG hr, far less than the viscous timescale. How- 



the disk edge, and fig 



^K{Rd) — ^sys is the beat 



ever, in the case Viallet & Hameury ( ^2007 ) studied they 
found that some of the gas may take longer, up to several 
system orbital periods, to reach the nozzle. 

3.2. Bifurcation (DCB) and warping estimates 

Matter streaming from the LI point typically collides 
with the edge of the disk, which usually sits outside the 
circularization radius due to viscous spreading. Because 
the disk is relat ively cold at this radius, the collision is 
highly ballistic ( Armitage & Livio |1998 ) . The fraction 
of matter streaniing around the disk instead of immedi- 
ately joining it can then be estimated simply by finding 
the altitudes at which the vertical disk and stream (both 
roughly Gaussia n) density profiles match, and supposing 
( Hessman|1999 1 that all the stream within this range im- 
mediately joins the disk while matter outside may stream 
further in. This yields a streaming fraction. 



fs{t) = erfc 



ln(fldo/pso) 

f - [Hs/HdY 



1/2 



(4) 



where p^o Pso are disk and stream densities at z = 
and the stream scale height Hg will not d iffer much 
from T 1 — we also refer to Hessman (19991 for fits to 



the results of Lubow & Shu (1975). 



Irradiation-driven warping of the o uter disk may also 
affect the streaming fraction. A gain, [Ogilvie fc Dubus| 



(2001) and Foulkes et al. (2010) suggest warping is pos- 
sible m LM(J and the latter work specifically finds 
a disk tilt of 10° likely for LMC X-3. However, both 
use an isotropic central luminosity, and the latter use an 
Eddington ratio in luminosity for LMC X-3 comparable 
to our derived maximum Eddington ratio in 9Jt^, so we 
consider their results an upper bound on warping. 

We generalize the fs (t) estimate to a stream that scans 
the edge of a disk tilted by an angle above the 

orbital plane. Here, the vertical density centroid follows 
zq = Rd sirL^d d(t) cos{ilsyat)) where Rd is the radius of 



frequency between the Keplerian frequency at the disk 
edge and the system orbital frequency. The finite travel 
time of the stream should add a roughly constant delay 
of order the local free-fall time, and for now we ignore 
this effect. Assuming d-dd/dt <^ ^syn, the altitudes of 
equal density are 



z± 



zoH,±Hdy/z^ + {H^ 



Hi) 



ln(p<io/Pso) 



(5) 



We will see that EC can depend very non- linearly on fs (t) 
at the boundary, but for now we use the orbit-averaged 
fs{t) as a gauge of plausible DCB strength: 



ifsit)) 



erf 





+ erfc 















(6) 



We plot {fs) at the outer boundary for relevant ranges 
of total supply rate, Mtot, and i?^, and for ^d of 0° 
and 10° in figjG] For an untilted disk, the contours are 
explained by the drop in disk scale height with radius and 
much slower drop with accretion rate, while for a tilted 
disk, the scanning greatly washes out the Rd dependence 
leaving accretion rate as the dominant factor. Our simple 
estimate also does not res olve the fate o f the s urviving 
stream beyond the edge (Foulkes et al. 2010 do, but 
unfortunately not for LMC X-3 in particular), but should 
bound the fraction of mass diverted. 
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Fig. 6. — Solid and dashed contours show {fs) for a disk edge 
tilted by 0° and 10° respectively, for a gas temperature of 16500K, 
and range of relevant Mtot E'nd outer disk radius. Tilting the disk 
edge generally increases (/s), but can also substantially change its 
dependence on the parameters, with the greatest effects at large 
Rd and Mtot- 



3.3. CHW Prescription and Estimates 
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In Begelman et al. 



(|1983|), the authors considered 
Compton 



an opticaUy thin corona subject to Compton heat- 
ing/cooUng (ignoring bremsstrahlung and other heat- 
ing/coohng mechanisms) and pointed out that accretion 
X-rays can heat the corona at all radii up to a tempera- 
ture, Tic, at which inverse-Compton heating and cooling 
equilibrate. Whether a wind is launched at a given ra- 
dius then depends mostly on whether this Tic is greater 
or smaller than the local virial temperature, Tvir, and the 
authors define a radius R^q by where the temperatures 
are equal, as well as a critical luminosity, L^j. -^Edd/"^^ 
at which the gas can be Compton-heated to the virial 
temperature within the sound-crossing time of the lo- 
cal corona's scale height. Because the tidal gravitational 
field falls off faster than the source luminosity, gas fiows 
out most easily at large radii. Though primarily a func- 
tion of source X-ray luminosity and radius, the shape of 
the source spectrum can af fect the mass-loss rate s lightly 
but we ignore this effect. 



Begelman et al. ( 1983 1 com 



puted mass-loss rates for total M^, < Mtot, wh ile later 



work addre sses dynamics and wind limit cycles (jShields 
et al.||1986|). 



Woods et ai. ( 1996 1 performed simulations to test the 
previous analytical prescription and amend it slightly — 
mostly by noting a shift in the location of R-^q and pro- 
viding corrections for low luminosities that do not im- 
mediately concern us. We take their fitting formula for 
wind losses per unit area 



dM, 



^ _■ 2/3 / 1 + (0.125r; + 0.00382)7C 
dA ""'"^'^ I 1 + (7^4(1 + 262^2))"' 



1/6 



(7) 



X exp 



-(l-(l + 0.25r^)"^''^)V2C 



where the normalization rhch is the ratio of corona pres- 
sure to sound speed at i?;^, ^ = R/R^q ~ 2R/R^-j.^, and 
their r] = L/L^^.. We then also introduce a factor f^h 
in r] = fxhL/L^j. for how well X-ray luminosity from an 
inner disk Compton heats the outer corona compared to 
the point source considered in the references. Although 
the outflow geometry may permit parts of the outflow 
to eventually reach low inclinations relative to the inner 
disk, the chief hurdle is heating the gas when it is sitting 
deepest in the tidal gravity field. Integrating cos i over 
the solid angle subtended by the outer corona versus half 
the disk's sky gives f^h ~ 0.025. This factor suppresses 
CHW considerably, while fxh ~ 1 implies CHW will have 
significant effects at the maximum observed luminosities 
(figlTl. Furthermore, the fxh for depleting disk flow is 
likely different and smaller than the corona as the X-rays 
will have to reach higher inclination, and heat conduction 
from a transition layer will be competing with advection 
by the wind. 

Winds may also be driven by other means, i.e., 
magneto- centrifugal an d line driving, but extensive simu- 
lations by Proga ( 2003 ) with parameters relevant to LMC 
X-3 indicate that these losses in LMC X-3 will be at most 
a few percent of the total accretion rate. 

To gauge CHW self-screening, or screening the com- 
panion, consider a wind carrying away lO^^g (total, 
half this per disk face) at the local sound speed at R^^. 
If the density did not fall off with radius, the Thomson 
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Fig. 7. — Prescription Mw losses per decade in R/R-^^ with left 
(right) vertical axis showing Eddington ratio when irradiation ef- 
iiciency /^.^ is 0.01 (1). The dashed (dotted) lines show a radial 
extent of 2R^^^^ and luminosity range for LMC X-3 with (without) 
reduced f^^- 



optical depth would be: 

0.5 X 10i9[g s-i]/mp 



ncTrpAs 



Rfcim^]kB/^p: 



1/2 



CTrnl a - 



Ric) ^ 0.15, 



(8) 

where a is orbital separation. That this extremely gener- 
ous upper bound gives marginal absorption indicates the 
Compton wind will not screen the companion. Instead, 
CHW and SRM will likely dampen each other's contribu- 
tion to Mtot -variability seen at inner radii as additional 
X-ray luminosity simultaneously increases supplied 
by the companion and lost to space. However, their 
interaction could enhance the scaling of M^/M^ with 
at large radii. 

3.4. PHII limits and discussion 

The PHII is fundamental to the picture of transient 
BHXRBs and thus to future extension of our work, but 
the physics itself is not trivial to implement let alone fully 
understood as the (60 page) review by Lasota (200l| 
attests. However, for LMC X-3, the strong, persistent 
disk emission should generally stabili ze th e disk within at 
least li^circ ^ill later show (^.1 1 that even dras- 

tic disk variability outside R^^j.^/25 is still too viscously 
dampened to explain observations, though the PHII may 
still contribute to the magnitude of disk variability, and 
likely plays an important role during complete state tran- 
sitions. 

Taking either the mean or median of Wlj^{t) over our 
data set as a suitable proxy for supply rate gives {M^} w 
O.lAfgjj^, and while the disk beyond ~ R^ij.^/3 will be 
cool enough to experienc e the PHII absent irrad iation 
at O.lMgjjj (e.g. fig.l of Janiuk & Czerny 2011), irra- 



diatio n can stabilize mo re and poss ibly aff 



ot the outer 



disk ( jCoriat et al.||2012| . Work by lbubus et al.| ( |1999 l 
indicates that LMC X-3 s disk would become susceptible 
to instability just beyond R^i^.^ at lO^^g s~^ for typical 
values of a (0.1) and an overall accretion to irradiation 
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efSciency factor C, originally fit to light curves of the 
BHXRB A0620-00 and roughly consistent with simple 
calculations based on an annulus-to-annulus irradiation 



et al. (2002) : 



geometry (see discuss ion in Kim et al. 1999 and com- 
parison at the end of |Uubus et al.| |1999| to [King et al.| 
[1997, ). 

Because we do not see 9jl^(<) decay on i?^.;^^ -viscous 
timescales in LMC X-3, it appears that the PHII would 
also lack a large span of starved inner disk for a heating 
front to propagate through. After the long, complete 
state transition of figjSlhowever, the disk recovery is flare- 
like, consistent with the notion that the PHII can play a 
significant role in LMC X-3 at low enough disk blackbody 
flux. 

3.5. EC Background and Implementation 

As stated in g] EC may occur if the disk and corona 
are thermally coupled — if the disk cannot efficiently radi- 
ate away corona heat conducted onto it, nor sufficiently 
cool the corona via inverse-Compton cooling, then it will 
experience net heating and evaporate, but otherwise it 
cools the corona which then condenses onto it. Thus 
the mass- exchange, or "EC" rate M^, is sensitively de- 
pendent on the balance of heating and cooling, and the 
very different scalings of heating and cooling mechanisms 
involved make possible a wide variety of behaviors. At 
present, several EC models incorporate viscous and com- 
pressive heating, bremsstrahlung, and inverse-Compton 
cooling in the accretion flow including LTMIIM07 and 
MP07. 

Besides separating the thresholds for disk for- 
mation/destruction normally degenerate under a 
bremsstrahlung-only density criterion via inverse- 



Compton coohng and thus engender hysteresis (Meyer 



Hofmeister et al. 2009), it is also possible to evaporate 
the outer disk but condense it back onto the inner disk 
rapidly enough to drive correlated rises (and falls) of 
3Jl^ with 9Jt^ (again, fig. 8 of MP07 and prominently 
in the left panel of our fig 11). It is also possible to 



preferentially evaporate the middle of a disk to t he point 
of destroying it as visible inlMeyer et all (|2007|), MP07, 



and several of our simulations. 

For our initial EC implementation, we do the follow- 
ing. We assume azimuthal symmetry for the accretion 
flow and divide it into 45 logarithmically-spaced radial 
zones with a single virtual corona zone associated with 
each disk zone (i.e. the code is 1.5D). We evolve the 
disk by solving mass fluxes with the standard viscous- 
disk equations (Frank et al. 2002j ) and a simple donor- 
cell scheme. Meanwhile we assume that the local corona 
properties and EC rates match those of the steady-state 
corona and solutions from LTMIIM07 for the same 
local accretion rate and evolve the corona working inward 
from the outer boundary condition. 

More specifically, for each disk zone we compute zone- 
boundary (j ±1/2) velocities 



3 {{vYsR^VtK)j - {vY.R^VtK)j+i) 



(9) 



with a viscosity based on a standard thermal equilibrium 
thin disk solution assuming Kramer's opacity as in , Frank, 



Vd = 2.13 X IQJ 



9 / ^^4h 



Mr- 



5/7 /QR\ 15/14 



The overall disk and corona evolution is then governed 



by: 



A(S],A^),/At=(Af,) 



(Aff- 



and 



(11) 



(12) 



where AA^ is the zone area. The fluxes are also limited 
so as not to draw mass from a disk zone or the corona 
flow than is physically available, and the "Hx" empha- 
sizes that we are plugging in the EC rates of LTMIIM07 
as a function of radius, and local and effective disk 
temperature. If and are anywhere comparable (of 
the same order of magnitude) to the viscous disk fluxes, 
then the computations are performed only once. Other- 
wise, the latter two steps are relaxed further, and allowed 
to change but not Af^, until their iterations converge 
within a given tolerance (or exceed an iteration limit). 
We note that large \MJ can lead to oscillatory behavior 
with this simple scheme, which can be subdued but not 
fundamentally fixed with smaller time steps and toler- 
ancea 
(fi 
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^hjs can be seen in the results of our simulations 
13 1 as small sudden jumps in Af^ (and somewhat 



in when EC is strong at small radii) , but by changing 
time step and tolerance we have found that this does not 
impact the general, longer timescale features that imme- 
diately concern us. We also explored small changes in the 
number of radial zones, obtaining similar results with 40 
or more grid points, but our results diverged very quickly 
for coarser grids. 

For steady input conditions, we confirm that our 
method does well reproducing cases considered by 
LTMIIM07 (treating the disk fully lets it spread viscously 
to larger sizes mildly enhancing condensation). To test 
our method's treatment of dynamic behavior, we looked 
at the same case MP07 studied: a disk spanning 3000 
Schwarzschild radii around a IOMq black hole with fixed 
boundary corona fraction fg of 0.1 and mass supply rate 
of 10"'^ times Eddington M-^^^. Their time-dependent 
method evolves the corona seff-consistently on its dy- 
namical timescale making it more physically realistic, 
but this also requires many more time steps. We found 
that with default parameters and physics our code never 
evaporates any part of the disk, while MP07 predict the 
formation of a gap that eats its way inward. However, we 
also found that if we scaled up the heat-conduction fluxes 
predicted by LTMHM07 for zones where inverse Comp- 
ton cooling does and does not set the electron temper- 
ature by a factor of three and flve respectively (adjust- 
ing the formula identifying the zones accordingly) then 
we do obtain good agreement with MP07 (see figjsj and 
their fig. 8). By preferentially scaling up heat fiuxes in 
the zones where inverse-Compton cooling limits electron 
temperature we obtained more rapid evaporation start- 
ing further inward while doing the same for heat fluxes 
in non-Compton zones led to increased stability. In this 
particular case, we saw little change when lowering the 
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magnetic-to-gas pressure ratio, /3c, a global constant in 
LTMHM07, from 0.8 to 0.1. 



or Compton-cooling (at high dJl^) limiting te mperatures 
may cross the coupling temperature (e.g. 
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Fig. 8. — The larger, darker empty squares and filled diamonds 
show OT^ and OT^ respectively at the inner boundary for a run with 
modified heat conduction fluxes (see text j ]3.5[ l, and the smaller, 
lighter symbols show the corresponding 9Jt , and QTt^ for default pa- 
rameters - here we plot time logarithmically for more direct com- 
parison with figure 8 of MP07. 

Both the models of LTMHM07 and MP07 necessar- 
ily neglect, or precede, some additional physical effects 
which may be relevant to our early results so we dis- 
cuss them here (and summarize i n fig [9| to motivate 
our more exploratory simulations (^4.2 and fig 13). In 
both models, condensation is a smooth, unresolve d flow, 
but applying the results of Wang et al. (2012) shows 
that the corona is liable to clump at radii greater than 
roughly lOORid under typical conditions for LMC X-3. 
Such clumping potentially increases cooling (thus con- 
densation) efficiency. Both LTMHM07 and MP07 use 
Spitzer electron conduction throughout the problem do- 
main, and both suspect that the effective thermal con- 
duction coefficient k may be significantly smaller. Al- 
though the degree of tangling in the magnetic fields of the 
transition zone is far harder to con strain, it is amenable 
to parameterization. Meanwhile, Cao (2011) provides 
a recent calculation for how much the ordered compo- 
nent of field shifts from predominantly poloidal outside 
~ lORid to predominantly toroidal inside. Lastly, nei- 
ther model includes mechanisms to spontaneously pro- 
duce corona, a poin t MP0 7 especially emphasize. Indeed, 
since 



Galeev et al. ( 1979| derived that within a certain 



radius, the buoyancy of magnetic loops formed within the 
disk can outpace their reconnection leading to a carpet 
of buoyant loops, this solar-like corona has often been 
invoked as a partial or complete so urce of corona. For 
LMC X-3, the condition on radius in Galeev et al. ( 1979 1 



gives R < 300(Afg(j^/9Jt^)i?s- It is hard to imagine this 
mechanism alone generating 971^ pulses that anticipate 
dJl^ drops lasting substantially longer than the viscous 
timescale at ~ IOO-R5, but it may play an important role 
by reheating and replenishing the corona, and affecting 
the local magnetic field geometry. 

An important caveat in applying the LTMHM07 model 
came to our attention after running our simulations. For 
> 0.1 Eddington near R ~100i?5, the conduction 



Bradley & 



Frank| 2009[) . Under these conditions, the model will tend 
to overpredict condensation and thus exaggerate the am- 
plitude and duration of the sympathetic mode, but trig- 
gering the effect requires strong, correlated 971^, to 
already be underway. 

4. SIMULATION RESULTS 
4.1. Models without EC 

To illustrate how a combination of non-EC mechanisms 
has difficulty explaining the magnitude of drops and 
especially the pattern of rapid decline versus slow recov- 
ery, we employ a very simple model where disk accretion 
is computed using the analytical machinery of ?. The 
latter is derived assuming a time-independent viscosity 
with power-law dependence on radi us, and we r e ferred 
to the thin-disk solutions on pg.93 of Frank et al. (2002 1 
for all relevant viscosity parameters. T'his method pre- 
vents incorporating radial dependence of M^, but be- 
cause wind losses fall rapidly with decreasing radius, and 
since we predict they are fairly weak anyway, assuming 
that they take place near the boundary does not inval- 
idate the main results of this toy model. This method 
als o pr ecludes incorporating he PHII, but as discussed 
in §3.4| the main effects of the PHII should typically be 
limited to radii beyond ^ ^circ LMC X-3. 

The left panel in fig[lO] shows a calculation with this 
reduced model geared toward reproducing the first disk 
drop in fig(5] We set fs manually, and SRM and wind 
losses were also computed beforehand as functions of 
the observed 9Jt^. To reproduce the depth of the first 
drop, we first allowed sustained of 100% and when 
this proved insufficient we moved the outer disk radius 
inward to a mere O.OAR^^j.^ for the runs shown, still em- 
ploying fs of 100%, and leading to massive OJtj, pulses 
compared to our estimates. This can be corrected some- 
what by invoking a wind stronger than our estimates. 

Besides requiring unrealistic values with respect to our 
estimates, and an extreme ad- hoc disk truncation, the 
toy model resists efforts to simultaneously improve agree- 
ment with other major features of the data. To improve 
model-data agreement for the second disk drop in the 
left panel of fig 10 without increasing disk radius (which 
would obviously undo agreement with the first 9Jl^ drop) 
requires either increased SRM or increasing the height 
and duration of the second coronal pulse. The latter 
will generate obvious disagreement with the second DJl^ 
pulse by attaching a tail that is very clearly not observed; 
the former significantly reduces the amplitude of the first 
coronal pulse relative to the second so that one must in- 
voke a stronger and more complicated wind mechanism. 

Again, the chief problem is that a disk flow will be vis- 
cously smeared too much to match observations unless 
variability is driven at a relatively small radius where 
even maximally efficient CHW should be insubstantial, 
and the PHII should not operate nor regularly drive heat- 
ing fronts. In the next section, we will show how EC 
may introduce a evaporation-to-condensation transition 
or gap at radii comparable to the outer edge of the arbi- 
trarily truncated disk of the toy model, but that it also 
tends to overpredict condensation at inner radii, gener- 
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Fig. 9. — A cartoon summarizing possible complications to the model especially concerning EC (j ]3.5[ j]5ll. At small radii, the magnetic 
field (dotted lines) in the corona and at the disk-corona interface may be significantly non-poloidal thus suppressing condensation (A), and 
further outward buoyant magnetic loops may still alter the magnetic field geometry besides introducing additional reconnection heating to 
the corona to continue supressing condensation (B). MHD winds might also be stronger than predicted and carry away more of the corona 
(C), while clumping of the corona at large radii may instead enhance condensation over evaporation there (D). The potential for the corona 
to viscously outflow radially wherever it achieves large density gradients (E) may also significantly affect our results. 
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Fig. 10. — The solid curves in the bottom (top) panels show 9Jt^ (QJ^c) simulated with the simplified model discussed in [4.1 and empty 
circles show the observed 931^ (Q^tc)' where the simulation units first are chosen to match the observed and simulated initial OJl^, as the 
simple model's disk machinery has no direct dependence on absolute accretion rate. The data points in the top panels are then scaled so 
that the maximum observed equals the initial which is a much larger absolute scale than our estimates suggest. The right panels 
are for a run with greater variability in the total mass supply which is shown as a dotted curve in all panels. The dashed curves show the 
effective corona/disk inputs at the outer boundary so that the remaining difference between solid and dashed curves in top panels indicates 
the wind loss. 



ating correlated fH^-OT^ rises and falls inconsistent with 
observations. 



4.2. Models including EC 

Except where specifically noted otherwise, for these 
simulations we again use ad = 0.1, but a corona ac — 0.2, 
the standard Spitzer coefficient for electron thermal con- 
duction, a /3c = 0.8 for the LTMHM07 EC prescription, 
the observationally favored R^^^^ — 2.7 x lO^^cm, assume 
Rid is 3i?5, and set^/j = 0.03. We note here that the 
SRM results from §3.1| specifically correspond to a value 
of -^"off/^*!eff = 0-4 for = lO^Sg s-\ and to 6.5 
scaleheights between the nozzle and the stellar surface 
for T'f = 165G0K. 

For the EC-inclusive model, we first simulated a disk 
with standard density profile extending to the circular- 
ization radius reacting to a mild Gaussia n pulse, 
and show the results in the left panel of fig[ll] Two im- 
mediately remarkable features include the sympathetic 
rise and fall of 9Jt^ with 93t^, and the general saturation 

of STtg response when the outer disk is most intensely si- 
phoning onto the inner disk. This saturation physically 
arises from the steep transition between evaporation and 
condensation. In the mass exchange model of LTMHM07 
that we employ, the local evaporation/condensation rate 

M^, varies with the local corona accretion rate in terms 



of Eddington ratio, to^,, like 

~aml/^(l-6mf/'') (13) 

where a and h are functions of many other parameters, 
local variables, and radius itself. If these other variables 
vary weakly with radius, then a very dense corona at 
some radius will lead to efficient condensation slightly 
further inward, and subsequent changes in about zero 
will be driven by the weaker variations in the critical 
value of . 

The issue of sympathetic accretion prompted us to con- 
sider a scenario in which the outer disk mass most vulner- 
able to being siphoned has already been evaporated away, 
such that there is a gap in the outer disk. We first studied 
the effects of simply truncating the disk, and show an ex- 
ample with radius 3 x lO^^cm and default LTMHM07 EC 
parameters in the right panel of fig |ll[ Truncating the 
disk to this radius prevents triggering sympathetic ac- 
cretion while generating appreciable COt^j variability and 
diminishing, but not eliminating, EC's role in amplifying 
variability over the simulation domain. 

We next attempted to generate this gap self- 
consistently, while preserving the interior aspect of the 
accretion flow that works fairly well. To this end, we first 
subjected a full disk to a constant at the boundary. 
The run confirmed that a small seed corona can rapidly 
evaporate the outer disk, that this corona is immediately 
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condensed only slight inward, but also that the process 
of forming a complete gap would take on the order of 
years in our standard EC implementation. To pursue 
this idea further, we ran simulations with a gap already 
inserted of which fig fT2] is representative. Besides the 
initial relaxation of DJl^ due to viscous spreading of the 
inner disk exceeding and resupply via condensation, one 
can see that sympathetic accretion is still an issue, and 
variability inward of the gap is severely suppressed 
again, as the high corona fraction of the gap triggers the 
saturation effect described above. 

Since our standard model and implementation of EC 
faces fundamental problems in reproducing major fea- 
tures of the data, we studied the effects of introducing 
physically motivated, if not yet rigorously justified mod- 
ifications. Thus far, it appears that the most successful 
modifications follow a fairly strict pattern of effectively 
raising the coronal heating and critical evaporation-to- 
condensation over the innermost two decades in 
R/ Rid- The latter is nearly inversely proportional to b in 
eqnJTS] which in the model of LTMHM07 scales with the 
corona viscosity parameter Oc, radius, electron thermal 
conduction coefficient k, and gas-to-total pressure ratio 
Pc as 

/ n \ -1/10 
6^/?e.V5^-l^A^(^^j (14) 

though we should point out that /3c enters their model 
strictly via a prescription for compressive heating. In 
the simulations producing fig |13[ we raised ac linearly 
from 0.2 to 0.4 inwards over lOO-Rid, and independently 
adjusted b over radial zones spanning 10°-10^'^. 10^'^- 



and lO^ '^-lO''-^ in R/R^d- For run A in fig 



13 



we 



10-^ 

scaled b in these zones by 0.1, 0.2, and 0.4 respectively, 
and for run B by 0.4, 0.6, and 0.8. Although the simu- 
lated pulses are far more massive than observations 
indicate, these modifications very clearly control the de- 
gree of hysteresis versus sympathetic accretion. If the 
overall magnitude of EC is smaller, then for higher fs, 
and/or SRM moderately larger than expected, this mod- 
ified model could reproduce the observed hysteresis. 

The most conceivably adjustable parameters in eqn 14 
are especially if understood to include other heat- 
ing mechanisms, and k. Since the current contrast be- 
tween simulation and observation still favors weaker EC, 
the requirement on additional heating would not likely 
be as extreme as implied by our modifications, but this 
then requires even greater deviation in k which affects b 
rather weakly, and these parameters are not necessarily 
independent in reality. 

5. DISCUSSION AND CONCLUSIONS 

Thus far, we still cannot offer a very definitive solu- 
tion for LMC X-3's behavior, but we have more rigor- 
ously examined several physical mechanisms popularly 
invoked to explain variability in LMC X-3 and addition- 
ally considered evaporation and condensation. We have 
found that if condensation is suppressed at inner radii 
(or over-predicted in the current model), then EC may 
naturally reconcile the observational evidence for RLO 
accretion and associated circularization radius, as well 
as the large amplitude, long duration, and rapid declines 
in hysteresis episodes with our estimates of other major 



viscous timescale variability mechanisms and the viscous 
dampening that they would undergo. 

However, we wish to emphasize that our current EC 
implementation and the steady-state theory informing it 
by default led to excessive condensation when compared 
as accurately as possible to observations, and it led to 
significantly different predictions for the p artic ular low- 
Mg case studied by MP07. As discussed in S 3.5 both the 
LTMIIM07 and MP07 models necessarily neglected some 
physics, some of which might enhance heating and/or 
suppress conduction closer to the black hole, and thus 
help explain the discrepancy between observations and 
the predictions of our code with the default EC prescrip- 
tion. 

To this end, we have reproduced and are testing a code 
that follows MP07 and evolves the disk and corona mass 
and energy equations self-consistently. An additional ad- 
vantage is that we can naturally include physics behind 
the PHII by incorporating detailed results for disk cool- 
ing as a function of density and central temperature from 
previous work - crucial to studying transient systems. 
However, this explicit method suffers from advancing by 
a very small time step. We have started building a par- 
allelized implicit method which we hope to develop fur- 
ther during simulations with the explicit method, but we 
also hope to find ways to save on excessive computation 
through better physical understanding of the problem. 

Our estimates for the reduced efficiency of CHW, and 
cursory examina tion of theory results for MHD winds in 
similar systems. Proga (20031 suggest that winds in gen- 
eral will have little attect on accretion dynamics in LMC 
X-3. However, we will continue to consider how efficiency 
of CHW might be increased, or that MHD winds may be 
stronger than anticipated (e.g. King et al. 2012|. 

Other assumptions in our current modeling that bear 
repeated mention include the alpha-prescription disk, 
and how we infer and interpret the disk and coronal 
accretion rates. Regarding the former, it is at least ex- 
pected that under conditions associated with jet flow, an- 
gular momentum transport via the magnetically driven 
outflow may become substantial or dominant compared 
to viscous transport (i.e. the magneto-rot a.tional insta- 
bility) at least within the i nner flow (e.g. Zanni et al. 
2007| |Casse fc Ferreira|2000l ). Evolution of the magnetic 
fields in the disk ma y also lead to intrinsic varia bility 
out to ^100i?i(i, as in'd e Guiran k, Ferreira| (2011 ). The 
potential for corona flow to stall centritugally as a func- 
tion of external flow conditions and local viscosity (e.g 



Chakrabarti fc Titarchuk|1995[|Garain et al.|2012| might 
be reafized in LMC X-3, but is also usuafly expected to 
occur well within the innermost IQQRid of the flow and 
to destroy the disk interior to the centrifugal shock, so 
it may be most relevant to state transitions. If more 
frequently prevalent though, the latter could alter our 
picture of spectral production, enhance SUt^, variability 
especially on shorter timescales, and change the dynam- 
ics of EC. Based on the pattern of the power-law com- 
ponent to anticipate declines in blackbody flux and the 
viscous recovery timescale of the blackbody component, 
we are still naturally inclined to favor a picture where 
variability is driven outside-in so that mechanisms like 
these and EC would predominantly accelerate and en- 
hance OJlj^ declines driven by known mechanisms operat- 
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Time [d] Time [d] 

Fig. 11. — Results with rough EC-implementation showing 9JT^ (empty black squares), OT^, (filled diamonds), Mtot (solid line), 
(dotted line), fs X lO^^g s~^ (dashed line) while without EC turned on (gray empty squares) is included for the right panel. The run on 
t he le ft uses the full circularization radius and shows strong condensation while the right panel shows a run with a disk size of 3 X lO^'^cm. 




Time [d] Radius [cm] 

Fig. 12. — Accretion history (left, see fig |ll| caption for symbol meanings) and evolution of the disk surface density profile (right) with 
snapshots at 80, 100, and 120 days shown m long-dashed, thick-dashed, and dot-dashed lines while the thin solid curves show the gap 
initial conditions relative to the standard thin disk. 




0""""2d""4d" 60 80 100 120 "l40 10' lO** 10'^ 10'" 10" 

Time [d] Radius [cm] 



Fig. 13. — The left panel shows accretion histories for OT^ (empty squares) and SK^, (diamonds) in cases A (thin, black), and B (thick, 

light gray) as described in j |4.2| where fs X 10^*g s~^ is also shown again, but and are o mitte d. The right panel shows evolution 
of disk surface density profile tor case A at 60, 90, and 120 days with the same convention as in fig |12| 



ing in the outer flow. 

Our immediate focus will be resolving the outstand- 
ing questions regarding the current mechanisms consid- 
ered, especially evaporation and condensation. After un- 



derstanding and constraining these better, we hope to 
extend our investigations to additional physics, systems, 
and phenomena, especially the transients and jet launch- 
ing. 
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APPENDIX 

RELATING ACCRETION RATES TO SPECTRA PRODUCTION 



For the blackbody spectrum, Zimmerman et al. ( |2005| note that XSPEC fits the maximum temperature and normal- 
ization constant to a temperature profile of the form T{R) = Ty^sj^{Rid/ R^^^- The fit to the peak temperature using 
this profile is only 5% smaller than the peak temp erature found fitting a temperature profile based on a zero torque 
boundary condition and color-correction factor f^^y ( Ebisuzaki et al.||l984[ ): 



Rid 

~R 



1/2^ 



1/4 



(Al) 



where 



1/4 



2.05T„ 



(A2) 



Because there is little difference in the fitted peak temperatures, because we plan to fix the black hole mass and Rid, 
and because we prefer the physically-motivated temperature profile we will use it instead. Integrating over disk annuli 
then gives the familiar formula for flux: 



F, 



mbb 



1 Airh COS i 



4 

col 



Rd 



Rid 



RdR 



exp(/iz///^„ifcBTeff(i?)) - 1 



(A3) 



Because Shimura & Takahara ( 1995 1 predict that /^.^j depends weakly on radius, accretion rate, and other parameters, 
we fix /^qj = 1.7. This then implies that the maximum disk accretion rate ranges from 0.07-0.29Mg^jj for Rid spanning 
1 to 3 Schwarzschild radii. We scale the IXft^ of figjs] by matching the observation with the highest blackbody fiux 
and temperature to the maximum 0.29M^^^. We note that if dissipation interior to the last stable circular orbit is 
significant, this will also put the actual 93t^ below our estimate (Beckwith et al. ]|2008 



Shafee et 



aI1|2008|. 



Infer ring 9Jl;,(t) requires additional assumptions but many are well constrained within ADAF theory (Narayan & Yi 
1995b). Specifically, theory predicts that corona ions are very effectively virialized at inner radii and much hotter than 
the electrons, whose exact temperature dep ends on many condition s, but is generally flat over the innermost lOORid and 
of order lOOkeV in the cases considered by Narayan & Yi ( 1995b). The former means that we can confldently predict 
scale height given corona ac while the latter provides some justification for choosing a constant electron temperature 
in corona emission c alculations. Taking R id = SRs, corona alpha parameter ac = 0.2, and gas-to-total pressure ratio 
Pc = 0.8, we obtain ( Narayan fc Yi|1995b ) the following estimates for corona density ric, corona height He, and coronal 
optical depth r in terms of mc = 3^c/-^Edd' 



nc(i?»d>c) « 1.2 X 10^^ mc{MQ/My^ J[cm' 



(A4) 



Hc{Rid,rnc)/Rid — 



(A5) 



r(i?j<i, mc) ~ 58mc 



(A6) 



However, we remind the reader that the latter is fairly sensitive to ac, scaling roughly like a^, ^. If we take the scattering 
fraction to be 7-",- = 1 — e~'^ then based on the scattering fractions returned by XSPEC for the mixed states, inversion 
gives us the simple SiH^ estimate shown as the solid line in figjsj 



We compare this simple estimate with a more detailed power-law flux calculation. Hua & Titarchuk ( 1995 ) find a 
Green's function for the output energy spectrum given the seed photon energy spectrum per scattered seed photon 
(specifically, their eqn.9), G,j{x,Xs,Te,pT)- The latter depends again on scattering fraction, as well as corona electron 
temperature T^, and the output (x) and seed (xs) dimensionless photon energie s (x^, = hu^, / kv,Tc) . Us ing the result 



of a more detailed calculation for the scattering fraction in a slab geometry from Zdziarski et al. (1994), 



1 



1 

27 



rEi(l,r), 



(A7) 
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we convolve G^{x,Xs,Tq^Vt) with the flux of seed photons that can scatter (overall Vt factor) out of the modified 
blackbody spectrum. In terms of r = R/Rid, and the power-law energy spectrum, Fe is given by: 

:(r3/4(l - 7-i/2)-i/4)/fc^T^) _ 1 (A8) 



47rd^ys Jo ./i /col exp{hiy, 



Note that we have assumed the corona emission is largely isotropic, but we have included the projection factor for 
blackbody emission from each annulus to half the height of the corona at r = 1, and we confirmed that r = 40 is 
a numerically acceptable cutoff. Fixing Tpn = 2.34 and = ISOkeV (Te dependence is relatively weak for Tg ^ 
max(T*,25 keV/fcg)) we tabulated integrated flux for a range of T^ax and 971^ to be inverted numerically, ultimately 
obtaining the points in the upper panel of flgjs] For relatively high T^^ and low t they agree fairly well with the 
simpler method, with the main difference due to the less step-like ^c~'T relationship in the more detailed Pr ■ However, 
to explain observations with higher Fpi and lower T^^,, the formula quickly requires excessively high 971^,, and at these 
implied optical depths the formula itself becomes unreliable. 
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